M-A 

k). J&fi-tSts 



, O 



LOW COMPRESSION RATIO, HIGH OUTPUT DIESEL 


(NAS A-CR— 169742) IGNITION AND COMBUSTION: N83-15309 

LOW COMPRESSION RATIO, HIGH OUTPUT DIESEL 
Final Report, Sep. 1979 - 30 Apr. 1979 

(Michigan Univ. ) 43 p HC A03/MF A0 1 Unclas 

1 G3/07 08039 


Final report of project conducted 
under NASA Grant No, NSG 3161 


Period of grant and report: September 1978 to 

April 30, 1979 


Work done by: 

/The University of Michig an 
College of Engineering 
Ann Arbor, Michigan 48109 

Project Directors: 

Professors Jay A, Bolt and William Mirsky 


July 1981 



- 2 - 


Table of Contents 

Page 

Mo. 

X. Introduction and Objectives , 3 

II. Summary of Results 3 

III. Equipment Changes and Improvements 11 

IV. Difficulties with and Changes made in Engine . . 14 

V. Computer Simulation of Diesel Cycle 16 

VI. Conclusions 26 

VII. Tables: 

1 - Best Performance Obtained at Cruise 

Power vs. Project Goals and Engine 

GTS 10- 5 2 OF 5 

2 - Comparison of Performance of 1977 and 1979 

Engine Configurations 10 

3 - Measured Engine Performance 24 

4 - Calculated Engine Performance 25 

VIII. Figures: 

1 - Engine Firing Time and Modifications, 

as Function of Date 6 

2 - Best Engine Performance of Experimental 

Single-Cylinder Engine 7 

3 - Best Engine Performance Under Earlier 

NASA Contract 8 

4 - Nozzle Hole Angles Related to Pistons 

AD-131 and AD-132 Top Contours 9 

5 - Ideal Diesel Cycle 18 

6 - Calculated versus Measured Indicated 

Specific Fuel Consumption 23 

7 - Calculated versus Measured Indicated 

Horsepower 23 

IX. Bibliography 29 

Appendix A - Title Sheet of Report No. CR-135300, 

Final Report for NASA Contract NAS-20051 

Appendix B - Computer Program Listings 



-3- 


I. Introduction and Objectives 

This project was undertaken as a sequel to an earlier 
NASA contract project entitled, "Lightweight, Low Compression 
Aircraft Diesel Engine," carriad out under NASA contract 
No. NAS 3-20051, under the direction of R. A. Kroegar of the 
Aerospace Department, College of Engineering, University of 
Michigan. Further details concerning this preceding project 
are given on the title sheet of the Report No. CR-135300, 
which is included herewith as Appendix A, and in reference 1*. 

This grant was made to determine the feasibility of 
converting a spark-ignition aircraft engine GTS 10-520 to 
compression ignition without increasing the peak combustion 
pressure of 1100 lbs/sq.in. The final engine contemplated 
was to utilize intake air heating at idle and light load 
and a compression ratio of about 10:1 with a small amount 
of "fumigation" {the addition of about 15% fuel into the 
combustion air before the cylinder) . 

The engine used was a modification of a Continental- 
Teledyne gasoline engine cylinder from the GTS10-520 supercharged 
aircraft angine. 

II. Summary of Results 

The normal process of development of the new test apparatus 
and of methods of test, coupled with some personnel problems, 
used the dollar budget ^ to such an extent that the combustion 
development remains about A months short of completion. 


*Numbers refer to those listed in bibliography. 

(1) The original cost estimate made at the University of Michigan 
to do the work was considerably greater than the amount NASA had 
available. Nevertheless, the work was undertaken on the basis 
that there would be a minimum of unplanned trouble. The test 
equpment "debugging", of course, was unexpected and this is 
explained on the basis that the cell needed to be revised more 
than expected because of increase in power for the 520 over the 
previous engine used in the test cell. 
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The reliable data taken on the last days of test 

indicate that the air utilization was poor, probably because 

10 ) 

the factors of air swirl' ' in the combustion chamber, 
injection of the fuel into the combustion volume, and fuel 
injection timing had not been optimized (note Figure 4) , page 9. 
See cruise power table comparison, Table 1, page 5 , which 
shows that 70% more air was required for equal power compared 
to the standard spark ignition engine whereas comparable 
well-developed diesels require about 33% more. We believe 
that a greater cut-and-try effort will yield the required 
improvement. The test conditions and performance for take-off 
power have not been evaluated in this work so far. The con- 
verted standard stark-ignition cylinder essentially withstood 
the testing without mechanical trouble even though some data 
points were obtained with 1300 lbs/sq.in. peak combustion 
pressure . 

As suggested above, there were many changes in equipment 
and improvements in the accuracy of the measurements as the 
work progressed. Some of these modifications are indicated on 
Figure 1, which shows hours of engine operating time as a 
function of calendar date. The result of all this is that 
we have little confidence in most of the data taken. For 
this reason only the performance data of Figure 2, and the 
tabulated comparison shown on Table 1, page 5 , have been 
included. These data were taken near the end of the grant 
period, and represent values which we believe to be most 
reliable. However, even these data points should be re-run 
to give confidence in them. 


(2) There is a considerable literature on the importance of 
air swirl and of air turbulence in utilizing air in the 
combustaon space and a recent reference is SAE 810477, especial- 
ly SAE 810255, and the list of references given in this paper. 
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Table 1 


Comparison of Continental Aircraft Engine GTSI0-520F 
and Best Performance of NASA Experimental Diesel Cylinder 

Cruise Power 


(75% x 435 = 324 BHP at 2900 

rpm = 170.2 

BMEP) 




GTS 10-5 20F 



Diesel 




Goal 

Present 

Compression ratio 

7.5 7.5 

7.5 

10 

9.62 

RPM 

2900 2900 

2900 

3040 

2600 

Fuel/Air 

.065^0.77 

.088 

.038 

.039 

Brake Sp. Fuel Cons. ,lbs/BHP-hr 

.43 .48 

.55 

,42 

.396(2) 

Indicated Fuel Cons. ,lbs/BHP~hr 



,39 

.356 

Peak Comb. Press, lbs/sq.in. 



1100, 

1200. 

Press. Rise Rate, psi/ca 



300 

186 

Brake MEP, lbs/sq.in. 



180 

188^ 

Indicated MEP, lbs/sq.in. 



200 

209 

Exh . Temp , °F 

1700 1600 

1500 

1500 

1355 

Intake Temp., °F 

192 192 

192 

200 

165°F 

Cyl. Head Temp., °F 



400 

386 

Smoke Hartridge 



10 

11.5 

Cooling Air Press. Drop "H 2 O 



6.0 

6.0 

Intake Air Press. In.Hg.Abs. 

36 36 

36 

53 

70 

Exh. Press. In.Hg.Abs. 



38 

50 

Air Utilization, lbs. Air /BHP. hr. 

C. 16 


8.8 

10.5 (2) 

Air Density, lbm/ft 2 




0.151 

Reference 

GTSIO-520-F 



July 


Report No. 629 


24-78 


14 July 70 





(1) Operation at this F/A is not recommended because detonation margin 
is small and because of higher cylinder head temperatures. 

(2) Mech. Eff. - 907 = - (based on friction HP curve for GTS10-520) 
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( 1 ) Scope, Needle Lift Trans., Poloroid, New Fuel Weigh , Replace Toledo with Link-Hartridge. 

Fig. 1. Engine firing time and Modifications, as function of date. 












INJ. TIMING SMOKE -H.S.U. IGN. DELAY COOLING AIR 




- 7 - 


Of POOR 


cn 

X 


CD 

X 

E 

u 


READINGS 169-180, log.sh. 28, 14 June 1978 
RPM-2600, IHP-51 (38KW) , IMEP- l79lbs/in 2 (1.23MPa) 
P M AX-t200 PSI (8,27 MPa), PRESS. RISE RATE-168 PSI/CA° 
INT. AIR PRESS. - 70“ Hg abs (177.8cm Hg (1.16 MPa /CA,) 



Fig. 2. Best performance obtained from GTS 10" 520 F 
modified diesel cylinder. 
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Fig. 3. Engine performance from Ref. I , earlier NASA contract. 


imer cylinder pressure 
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Contour of 




Fig. 4. Nozzle hole angles related to pistons AD-131 and AD-132 top contours. 


2.74 
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The best results at the cruise condition obtained 
from our single-cylinder tests are shown on Figure 2. This 
is in comparison with the data shown on Figure 3, which is 
a reproduction of Figure 29 , taken from reference 1, Report 
No. CR-135300. Table 2 below compares some of the significant 
values. It will be noted that substantial gains were accomplish 
ed in improving fuel economy and reducing smoke. 

The results of the computer simulation programs are given 
in a following section. 


Table 2 

Comparison of Single Cylinder Test Engine Performance 
at Maximum Output of 1977 Contract (Ref. 1) with 
Present Grant Best Performance Conditions: 


RPM 2600. 

Intake air temper ature,°F 225 

Intake air pressure 70" Hg. 

Exhaust pressure 50" Hg . 

Test run No. 174 


Ind. 

M.E.P. 

1977 1979 

175 179 


Ind. 

SFC (lbs/ind. HP-hr ) 
1977 1979 

.51 0.356 


Smoke Hartridge 
H.S.U. 

1977 1979 

80 11 . 5 
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III. Equipment Changes and Improvements 

During the period of this grant many test cell modifica- 
tions, equipment changes and improvements have been made. 

See Figure 1 - Engine Firing Time vs. Date (for the date when 
these changes were made) . The list of major changes and 
modifications which follow are indicative of the problems 
which occurred. 

1. The pressure transducers for the cylinder pressure 
(Kistler-type 6005) and for the fuel line pressure 
(Kistler-type 607 FX) were recalibrated several 
times. 

2. The two- trace oscilloscope was replaced by a four- 

trace Tektronix oscilloscope to allow simultaneous 
recording of: cylinder pressure, fuel line pressure, 

needle lift, and crank angle degree marks. 

3. The needle lift instrument installation was improved 
following a visit to the Stanadyne Co., Hartford, 
Connecticut, where many details were discussed. The 
previous installation used during the earlier contract 
resulted in broken injection springs and plugging of 
several nozzle holes of the injector AD-129 which was 
in the engine at the start of this grant. The broken 
spring is suspected to be the reason for the extreme- 
ly high smoke (70 and 80 H.S.U) and high specific 
fuel consumption (0.51 and 0.52 lb/hp-hr) of Figure 3 
from Report CR-135300) . 

4. The appropriate Polaroid camera was reconditioned 
and installed on the oscilloscope. 

5. It was found that the scale beam readings were more 
steady without use of the engine automatic speed 
control device. A spring was installed in the pump 
rack control linkage to insure constant rack 
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position during test runs. To measure rpm, a total 
revolution counter was synchronized with the fuel 
weighing timer. 

6. A new fuel weighing device was incorporated into 
the equipment set-up to improve the precision of 
fuel measurement. 

7. The leak-off of the injector was ducted back to 
the fuel measuring beaker, to avoid error from 
such leakage. 

8. Two additional temperature indicators were installed 
for direct measurement of cylinder head and exhaust 
ga-j temperature. 

h. A new air pressure regulating valve was installed 
in the air supply system which resulted in less 
fluctuation of the engine intake air pressure. 

10. The cylinder pressure transducer was mounted flush 
with the inner combustion chamber wall by means 

of a new adapter. 

11. The thermocouple probe for the exhaust temperature 
was found to be burned and was replaced. 

12. The Toledo gear-weighted pointer dynamometer load 
measuring unit was replaced by a Link-Unibeam 
System to (a) eliminate need for observer to average 
pointer reading variation due to vibration of point- 
er, (b) increase scale length when reading engine 
friction, (c) to prevent destruction of the gear teeth 
of the Toledo unit. 

13. A fuel leak near the base of the fuel line pressure 
transducer holder was repaired. 
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14 . The Hartridge smokemeter showed inconsistent 
readings and was reconditioned by the factory 
authorized service. 

15. A new AVL water-cooled quartz pressure transducer 
No. 12QP505 with special temperature compensation, 
and with a shield, was installed for reduced 
radiant heat effects. The pressure unit was also 
suspect because of zero pressure uncertainty. 

16. At higher engine speeds, the flexible exhaust hose 
failed several times and it was necessary to 
modify the exhaust system by lowering the exhaust 
surge tank and thus eliminate the bending and 
resultant stress of the flexible hose between the 
cylinder and the tank. 
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IV, Difficulties with and Changes made in Engine 

During the period of November 1, 1977 through July 31, 
1978, several engine-related difficulties were encountered, 
and corresponding changes and/or improvements were made, as 
follows : 

1. After several introductory runs of the engine in 

search of the cause of the very dark exhaust ( 70 

Hartridge units) and poor performance, injector 
AD-129 was taken out and dismantled. The injector 
spring was broken in four pieces and several spray 
holes were plugged with carbon. The manufacturer 
of the nozzle, Stanadyne, suggested that the most 
probable cause of that failure was improper installa- 
tion of the needle lift instrumentation. The nozzle 
was not repairable due to long operation with the 
broken spring. 

2. Mechanical pointers with positive locking devices 
for the injection timing and rack 1'nkages were 
fabricated and installed on the Bosch APE1B fuel 
pump, which eliminated slight creeping of the 
hydraulic controls. 

3. The high pressure fuel line (Bosch pump-injector) 
was replaced with one having 0.078 in. inside 
diameter instead of previous one with 0.083 in. 
inside diameter. This reduced some minor after- 
injection. 

4. When rn effort was made to increase the load of the 
engine at cruise power output, after-injection was 
again encountered. After experimenting with four 
delivery valves we found that the valve (Bosch 
special desigr X2708) with retraction volume of 

65 mm^ eliminated the after injection and all test 
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runs after mid-April of 1978 were done with that 
valve. All prior test runs were done with the 
Bosch valve X2615 (50 mm^ retraction volume) . 

5. After a short period of engine running at cruise 
power output, excessive blowby was observed. 

Teardown of the engine revealed that piston AD-131 
was badly scuffed. Teledyne-Continental Motors 

in Mobile, Ala. examined the piston-cylinder assembly, 
and found that the top piston ring had only 0,001 in. 
side clearance instead of the normal clearance of 
0.004 in., and concluded that this was the main 
cause of the piston scuffing. It was also suggested 
that we reduce the upper piston land diameter of 
5.216/5 I' 1 5 inches by about 0.005 inches, due to the 
piston being slightly higher than the original 
GTS10-52G piston and possible higher temperatures at 
the piston top. 

6. The support or retainer of the rocker arm shaft on 
the exhaust valve failed four times. After the 
first failure it is possible that the Helicoil was 
damaged and that all subsequent failures occurred 
because of the damaged Helicoil thread which caused 
low bolt holding force and resulted in a high stress 
range in the bolt threads. Teledyne-Continental 
Motors does not have a history of such failures, and 
it is believed that the fault should not be ascribed 
to the diesel combustion. 
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V. COMPUTER SIMULATION OF DIESEL CYCLE 

A sub-task of this study was to develop a simple 
computer simulation of the diesel cycle with the objective 
that it would aid in correlating the experimental data and 
could then be used to extrapolate to conditions which indicate 
improved engine performance. 

Two programs were developed during this study. The 
first, DIESELA, gives predicted engine performance parameters 

IMEP indicated mean effective pressure, 

ISFC indicated specific fuel consumption, 

IEFF indicated thermal efficiency, and 

IHP indicated horsepower 

based on input values for 

A/F air-fuel ratio, 

CR compression ratio, 

Tl temperature at the start of compression, 

PI intake manifold pressure, and 

P6 exhaust pressure. 

A limiting peak pressure of 1200 psia and a compressic ratio 
of 9.62, which correspond to the experimental test conditions, 
were used in all cases examined. 

A second program, DIESELB, was started in order to 
improve the simulation of the constant volume and constant 
pressure combustion processes. It is based on calculated 
equilibrium conditions for the ends of the ideal constant 
volume and constant pressure combustion processes. The 
objective was to calculate the thermodynamic properties at 
these states rather than rely on the use of average specific 
heats derived from combustion chart data. 
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1. Program DIESELA 

This simulation is based on a semi-empirical model which 
starts with the ideal cycle shown in Fig. 5 and then applies 
semi-empirical factors in an attempt to bring the model 
results into better agreement with measured engine data. 
Details of the model are as follows: 

(a) The temperature at the beginning of the compression 
stroke (Tl) is assumed to be equal to that measured in 
the intake manifold. The corresponding pressure (Pi) is 
taken to be equal to the intake manifold pressure (PI) 
times the volumetric efficiency. The initial charge is 
assumed to consist only of air. 

(b) Mass of the initial charge at state 1 is obtained from 
the ideal gas equation and the assumed pressure and 
temperature conditions from (a) above. The volumetric 
efficiency equation is taken from Taylor and Taylor 
(vol. 1, p. 157) . This equation was tested using 
calculated and measured air mass flow rates and the 
resulting equation shows the influence of engine RPM. 
Attempts to find effects of other parameters on 
volumetric efficiency showed no clear correlations. 

(c) Compression (1-2) is assumed to follow the poly tropic 
process PV n = constant, where the value for the 
polytropic exponent was obtained from Continental 
Motors Report No. 1069, Aug. 1967. The report contains 
a plot which correlates data from a large number of 
engines and gives values of the polytropic exponent as 
a function of both maximum pressure and compression 
ratio . 

(d) Heat addition at constant volume is calculated using: 


Qv = M*Cv (aver) * (T3 - T2) . 
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Pig. 5. Ideal Diesel Cycle 
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An adiabtic process is assumed and Cv(aver) is computed 
from combustion chart data for the temperature range T2 
to T3, the value of T3 being dependent upon the 
selected maximum cycle pressure P3. The amount of mass 
(M) is that of the air charge plus the amount of fuel 
required to reach maximum pressure during the constant 
volume process. An iterative process is used to get 
compatible values for T3 and this fuel mass. 

(e) Heat addition at constant pressure is used to determine 
state (4) using: 

Qp = M*Cp (aver) * (T4 - T3) , 

where Qp is equal to the difference between the 
available heating value of the fuel and that released 
during the constant volume process. An average value 
for Cp is computed from combustion chart data for 
the temperature range T3 to T4, using an iterative 
procedure similar to that used in (d) above. 

{f) Expansion (4-5) is assumed to follow a polytropic 
process having the same exponent used for the 
compression process . 

(g) The remaining processes are: 

(5-6) constant volume expansion to the pressure P6 in 
the exhaust manifold. 

( ) constant pressure exhaust process. 

( ) constant pressure intake process. 

(h) Indicated work for the ideal cycle is compared with 
the indicated work obtained from experimental brake and 
friction work values. The ratio of the two indicated 
work quantities is defined as the "diagram factor". 
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starts at the beginning of the constant volume combustion 
process using results from DIESELA. These include: 

T2 temperature at the end of compression. 

P2 pressure at the end of compression. 

V2 cylinder volume at the end of compression. 

NA moles of air in the fresh mixture. 

The program then computes the state of the combustion 
gases resulting from adiabatic constant volume combustion and 
assumed equilibrium at the limiting peak pressure (P3). The 
fuel required for this process (FV) is calculted and the 
remaining fuel (F-PV) is then consumed during the constant 
pressure combustion process. Equilibrium at the end of the 
constant pressure process is also computed. 

The program listing for DIESELS, the corresponding list 
of symbols and a sample output are shown in the Appendix B. The 
sample output gives the following information: 

CR compression ratio. 

T2 temperature at the end of compression in deg. K. 

A/F air-fuel ratio. 

PMAX limiting maximum pressure, (1200 psia) . 

For the constant volume process: 

Y(j) moles of specie J at equilibrium. 

X ( J) mole fraction of specie J at equilibrium. 

YTOT total moles of products at state 3. 

FUEL moles of fuel used during the constant volume 
combustion process. 

For the constant pressure process: 

Y { J) moles of specie J at equilibrium. 

X ( J) mole fraction of specie J at equilibrium. 

YTOT total moles of produsts at state 4. 
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Diagram factors for forty-eight engine runs were 
obtained and correlations with various operating 
variables were tried. Only mixture ratio was found to 
show a reasonable correlation. 

(i) Application of the diagram factor to the ideal cycle 

calculations gives predicted performance of the engine 
under various operating conditions . 


The program listing of DIESELA is given in the Appendix 
together with a listing of the nomenclature. The sample 
output/ shown in Table IV, gives values for the following 
input engine conditions : 

Pmax 
RPM 
A/F 
CR 

Computed values listed are: 

T5 Deg F. 

MV/MF ratio of mass of fuel used during constant 
volume to the total mass of fuel per cycle. 

DF diagram factor = {calc IHP/exper XHP) . 

IMEP indicated mean effective pressure. 

ISFC indicated specific fuel consumption. 

IEFF indicated thermal efficiency. 

I HP indicated horsepower. 


Tl 
PI 
P 6 


2. Program DIESELB 

This program attempts to improve the cycle simulation of 
the constant volume and constant pressure combustion processes 
by eliminating the need for combustion chart data. The program 
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VOL gas volume at state 4. 

T4 temperature at state 4 in deg K. 

PRESS pressure at state 4 in atmospheres. 

Some additional work would be required to incorporate 
this program in the DIESELA simulation. 

3. Computer Study Results 

Results showing calculated versus measured XSFC are 
given in Fig. 6 and similar results for IHP are given in Fig. 
7. Forty-eight test runs, taken at the end of the test 
program, were used in this comparison. Measured values are 
given in Table III while the corresponding computer results 
are shown in Table IV. The ranges of engine parameters covered 
in these test runs are: 


CR 

9.62 

(Fixed) 

RPM 

2405 

- 3007 

A/F 

19.15 

- 32.90 

Tl 

142 - 

286 F 

PI 

59.4 

- 79.4 " 

P6 

44.4 

-49.4 " 



calculated. ISFC 


-23- 


ORIGSMAL PAGE f9 
OF POOR QUALITY 



measured-lSFC 


Fig. 6 - Calculated versus Measured Indicated 
Specific Fuel Consumption 



Fig. 7 - Calculated versus Measured 
Indicated Horsepower 
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Table 3 - Measured Engine Performance 


RUN RPN fl/F Tl PI P6 INEP ISFC VE IHP 

169 2606 26.32 204 69.60 40.60 176.3 0.363 0.957 50.24 

170 2600 26.30 204 69.60 40.60 176.5 0.363 0.950 50.16 

171 2600 27,74 252 69.60 48.60 170.0 0.359 0.961 50.60 

172 2600 27.82 252 69.60 48.60 177.8 0.358 0.962 50.53 

173 2596 29.27 225 69.60 40.60 179.4 0.356 0.976 50.92 

174 2598 29.33 226 69.60 40.60 179.1 0,356 0.977 50.87 

175 2601 30.74 195 69.60 48.60 170.6 0.354 0.970 50.79 

176 2602 30.59 195 69.60 48.60 180.1 0.353 0.970 51.24 

177 2606 31.74 165 69.60 40.60 181.6 0.349 0.957 51.73 

178 2605 31.77 164 69.60 48.60 181.6 0.348 0.956 51.72 

179 2598 32.90 142 69.60 48.60 101.5 0.350 0.959 51.56 

180 2598 32.90 142 69.60 48.60 181.5 0.350 0.959 51.56 

190 2608 31.50 164 69.20 49.20 181,5 0.351 0.959 51.76 

191 2604 31.73 164 69.20 49.20 181.0 0.350 0,960 51.51 

192 2008 32.59 165 69.20 49,20 172.0 0.348 0.935 52.00 

193 2804 32.89 162 69.20 49.20 172.1 0.345 0,931 52.76 

194 2405 30,92 165 69.20 49.20 100.5 0.351 0.980 49.57 

195 2409 30.96 165 69.20 49.20 189.1 0.351 0.984 49.80 

196 2610 21.96 166 69.30 49.30 221.2 0.409 0,952 63.12 

197 2607 22.02 166 69.30 49.30 221.1 0.409 0,953 63.00 

198 2809 24.00 166 69.30 49.30 202.9 0.384 0.924 62.31 

199 2011 24.67 lt>6 69.30 49.30 203.1 0.383 0.919 62.41 

200 3005 25.37 165 69.30 49.30 191.8 0.300 0.885 63.02 

201 3007 24.98 165 69.30 49.30 191.7 0.386 0.-983 63.00 

200 3005 25.39 165 69.40 49.40 191.5 0.381 0.885 62.91 

201 3007 25.03 165 69,40 49.40 191.7 0,386 0.8B3 63.00 

202 2010 25.39 165 69,40 49.40 197.4 0.300 0.909 60.65 

203 2805 25.52 165 69.40 49,40 197.8 0.378 0.910 60.65 

204 2599 25.45 166 69.40 49.40 209.0 0.374 0.952 59.39 

205 2603 25.29 166 69.40 49.40 208.7 0.374 0.944 59.39 

206 2601 27.94 166 79.40 47.90 211.7 0.375 0,928 60.18 

207 2601 27.95 166 79.40 47.90 211.7 0.375 0.927 60.19 

208 2602 23.72 166 69.40 44.40 211.2 0.393 0,942 60.09 

209 2605 23.78 165 69.40 44.40 211.6 0.391 0.939 60.27 

210 2599 19.15 166 59.40 44.40 211.3 0.414 0.936 60.05 

211 2601 19.18 166 59.40 44.40 211.0 0.414 0.935 60.01 

212 2405 26.71 165 00.50 49.50 225.7 0,309 0.920 59.35 

213 2412 26.72 165 00.50 49.50 226.3 0.387 0.917 59.68 

214 2605 23.27 165 67,50 44.50 211.9 0.392 0,950 60.36 

215 2610 23.24 164 67.50 44.50 211.8 0.392 0.947 60.44 

216 2590 22.85 195 69.00 44.50 210.1 0.394 0.953 59.67 

217 2602 22.85 195 69.00 44.50 210.0 0.394 0.951 59.75 

218 2600 23.19 224 72.00 44.50 210,5 0.393 0.967 59.04 

219 2603 23.17 224 72.00 44.50 210.3 0.393 0.965 59.87 

220 2600 23.35 255 75.20 44.50 211.6 0.390 0.972 60.15 

221 2603 23.37 255 75.20 44.50 211.5 0.390 0.972 60.19 

222 2602 23.43 286 78.40 44.50 212.2 0.390 0.978 60.36 

223 2603 23.49 285 78.40 44.50 212.4 0.388 0.974 60.44 



Table 4 

P3= 1200 
RUN RPfi fl/F 


169 2604 26,32 

170 2600 26.38 

171 2600 27.74 

172 2600 27.82 

173 2596 29.27 

174 2598 29.33 

175 2601 30.74 

176 2602 30.59 

177 2406 31.74 

178 2605 31.77 

179 2598 32.90 

180 2590 32,90 

190 2608 31.50 

191 2604 31.73 

192 2808 32.59 

193 2804 32.89 

194 2405 30.92 

195 2409 30.96 

196 2610 23.96 

197 2607 22.02 

198 2809 24.80 

199 2811 24.67 

200 3005 25.37 

201 3007 24.98 

200 3005 25.39 

201 3007 25.03 

202 2810 25.39 

203 2805 25.52 

204 2599 25.45 

205 2603 25.28 

206 2601 27.94 

207 2601 27.95 

208 2602 23.72 

209 2605 23.78 

210 2599 19.15 

211 2601 19.18 

212 2405 26.71 

213 2412 26.72 

214 2605 23.27 

215 2610 23.24 

216 2598 22.85 

217 2602 22.85 

218 2600 23.19 

219 2603 23.17 

220 2600 23.35 

221 2603 23.37 

222 2602 23.43 

223 2603 23.49 
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Caleulated Engine Performance 


CR II T5 PI P6 HV/HF DP IHEP ISFC JEFF IHP 


9.62 284 1571 69.60 48.60 0.42 0.821 181.2 0.377 37.08 51.63 

9.62 284 1569 69.40 48.60 0.42 0.821 181.2 0.377 37.12 51.50 

9.62 252 14BB 69.60 48.60 0.42 0.833 184.1 0.369 37.95 52.33 

9.62 252 1434 69.60 48.60 0.42 0.834 1B3.B 0.368 37.99 52.24 

9.62 225 1408 69.60 48.60 0.42 0.847 186.0 0.359 38.90 52.80 

9.62 226 1406 69.60 48.60 0.43 0.847 i85.5 0.359 38.93 52.68 

9.62 195 1333 69.60 48.60 0.42 0.861 189.6 0.351 39,84 53.90 

9.62 195 1339 69.60 48.60 0.42 0.859 190.0 0.352 39.74 54.05 

9.6 1 165 1279 69.60 48.60 0.41 0.871 195.4 0.345 40.50 55.67 

9.62 164 1278 69.60 48.60 0.41 0.871 195.7 0.345 40.52 55.73 

9.62 142 1126 69.60 40.60 0.41 0.884 199.8 0.339 41.29 56.76 

9.62 142 1226 69.60 48.60 0.41 0.884 I99.fi 0.339 41.29 56.76 

9.62 164 1285 69.20 49.20 0.42 0.869 194.7 0.347 40.28 55.52 

9.62 164 1277 69.20 49.20 0.42 0.871 194.2 0.346 40.43 55.29 

9.62 165 123B 69.20 49.20 0.47 0.880 184.6 0.342 40.94 56.68 

9.62 162 1227 69.20 49.20 0.47 0.884 1B4.9 0.340 41.15 56.67 

9.62 165 1318 69.20 49.20 0.37 0.863 204.0 0.350 39.92 53.63 

9.62 165 1317 69.20 49.20 0.37 0.863 203.7 0.350 39.95 53.64 

9.62 166 1782 69.30 49.30 0.29 0.787 237,2 0,408 34,20 67.69 

9.62 166 1778 69.30 49.30 0.29 O.7B0 237.0 0.407 34.32 67.54 

9.62 166 1583 69.30 49.30 0.36 0.809 213.7 0.388 36.07 65.63 

9.62 166 1590 69.30 49.30 0.36 0.808 214.3 0.388 36.00 65.86 

9.62 165 1539 69.30 49.30 0.39 0.813 204.9 0.384 36.42 67.31 

9.62 165 1561 69.30 49.30 0.39 0.810 206.7 0.386 36.19 67.96 

9.62 165 1538 69.40 49.40 0.39 0.813 205.0 0.384 36.42 67,36 

9.62 165 1558 69.40 49.40 0.38 0.811 206.7 0.386 36.21 67.96 

9.62 165 1549 69.40 49.40 0.36 0.813 211.3 0.384 36.42 64.92 

9.62 165 1542 69.40 49.40 0.36 0.814 210.9 0.383 36.50 64.67 

9.62 166 1560 69.40 49.40 0.34 0.814 218.3 0.384 36.43 62.03 

9.62 166 1569 69.40 49.40 0.33 0.812 219.0 0.385 36.33 62.32 

9.62 166 1483 79.40 47.90 0.25 0.835 238.4 0.369 37.87 67.80 

9.62 166 1482 79.40 47.90 0.25 0.835 238.4 0.369 37.87 67.79 

9.62 166 1665 69.40 44.40 0.31 0.800 230.5 0,392 35.67 65.57 

9.62 165 1660 69.40 44,40 0.31 0.801 230.4 0.392 35.70 65.61 

9.62 166 1943 59.40 44.40 0.36 0.767 224.3 0.424 33.01 63.74 

9.62 166 1940 59.40 44.40 0.36 0.767 224.0 0.423 33.03 63.70 

9.62 165 1564 80.50 49.50 0.20 0.824 254.9 0.380 36.77 67.03 

9.62 165 1563 80.50 49.50 0.20 O.B24 254.6 0.380 36.78 67.13 

9.62 165 1681 67.50 44.50 0.32 0.797 226.8 0.395 35.43 64.61 

9.62 164 1683 67.50 44.50 0.32 0.797 227.2 0.395 35.41 64.83 

9.62 195 1730 69.00 44.50 0.32 0.794 224.4 0.397 35.19 63.74 

9.62 195 1730 69.00 44.50 0.32 0.794 224.2 0.397 35.19 63.79 

9.62 224 1737 72.00 44.50 0.30 0.796 222.7 0.395 35.40 63.29 

9.62 224 1738 72.00 44.50 0.30 0.796 222.7 0.395 35.39 63.37 

9.62 255 1759 75.20 44.50 0.28 0.79B 222.0 0.394 35.50 63.11 

9.62 255 1757 75.20 44.50 0.29 0.798 221.0 0.394 35.52 63.12 

9.62 286 1786 7B.40 44.50 0.27 0.798 221.7 0.393 35.55 63.08 

9.62 2B5 1781 78.40 44.50 0.27 0.799 221.7 0.393 35.59 63.08 


ORIGINAL PAG£ IS 
OF POOR QUALITY 
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VI. Conclusions 

1. Definite conclusions concerning the relative 
merits of the diesel cycle with low compression 
ratio and heated air for aircraft applications 
cannot be drawn from this work. 

2. Improvement in the air utilization of the 
experimental diesel cylinder must be accomplished 
before the merits of the diesel cycle can be 
evaluated in comparison with the present Otto 
cycle aircraft engines. 

3. The computer programs developed give results which 
are in general agreement with the experimental 
remits. The computer values of indicated specific 
fuel consumption are somewhat better than the 
computer values of indicated horsepower, the latter 
being high by about six percent. 
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VII, Tables 

Table 1 - Comparison of Performance at Cruise Power 
of Best Performance Obtained, vs. Goal, 
vs. GTS 10- 5 2 OP 

Table 2 - Comparison of Performance of 1977 and 1979 
Engine Configurations 

Table 3 - Measured Engine Performance 

Table 4 - Calculated Engine Performance 
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Fig. 1 - Engine Firing Time and Modifications, 
as Function of Date 

Fig, 2 - Best Engine Performance of Experimental 
Single-Cylinder Engine 

Fig. 3 - Best Engine Performance Under Earlier 
NASA Contract 

Fig. 4 - Nozzle Hole Angles Related to Pistons 
AD-131 and AD-132 Top Contours 

Fig. 5 - Ideal Diesel Cycle 

Fig. 6 - Calculated versus Measured Indicated 
Specific Fuel Consumption 

Fig. 7 - Calculated versus Measured Indicated 
Horsepower 
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APPENDIX B 

Computer Program Listings 



10 REH ‘DIESELA" 

20 REH DIESEL CYCLE ANALYSIS 
30 REH CV AND CP EONS DERIVED FROM COMBUSTION CHART 
40 DEFINT G.I-K 
45 LPRINT CHR4I15) 

50 R=53.34:P3 a l2O0U44:Cfl=9.62 

55 PI=P3/144000 

40 READ RN r RPH,Tl,PI,P6,Al 

80 P1=PI 170. 7264 

90 P6=P6»70.7264 

120 T1=T 1+460 

170 El=2,*912-8.6669E-4IRPH+2. 1226E-74RPHI2-1.9B59E-11IRPKC3 
190 Xl c (12563+1089tPX-8tPX[2-171IPX[3+35*PXI4)/lE4 
200 X2 a (17+9BQ6IPX-B783IPXt2+3113IPXI3-3B7tPXI4)/lE6 
240 IFIT>0THEN320 
270 LPRINT"P3=*jP3/144 

290 LPRINT 'RUN RPH A/F CR T1 T5 PI P6 hV/MF 
310 F0RI=lTO7fi:LPRINT“-“; sNEXT : LPRINT" 1 
320 PRINT’COHPUTING'iPRINT" ■ 

340 VE=E1I(.2B5+<CR-(P6/PI)]/(1.4I(CR-1)J) 

350 P1=VEIPI 

360 Vl a .050lltCRnCR-l) 

370 X=Xl-X2t ICR-6) 

390 AF=AI 

400 DF=.58O58+.O14B17»AF-3.9079E-4IAFI2+6.9442E-6*AF[3 
410 IT=IT+1 
420 FA=1/AF 

430 HA=P1IV1/ IRIT1) :HF=FAIHA:H=HA+HF 
440 B=«FI 18200 
450 P4=P3:V5=V1 

470 Hl=Pimi(CR[(X-l)-l)/(l-X) 

480 V2=V1/CR:V3-V2 
500 P2=P1ICRIX 
510 T2=T1ICRE (X-l) 

520 FR[l)=.2iFRI2)=.3 
522 F0RI=1T02 

525 HV < I ) =FR I IHHF : T3-P3IV3/ <RI IHA+HV (I) ) ) 

530 T=(T2+T3)/2 

540 CV=. 17276+4. 4512E-5IT-B.5504E-9ITE2+1.1852E-12ITI3 
542 ER i I)= (HV 1 1 ) +HA) tCVt (T3-T2) -HV 1 1 >118200: NEXT 

544 IFA8StER(2)KlE-5THEN560 

545 FR=(ERimFR(2)-ER(2)»FR(l))/(ERU)-ER(2)) 

547 FR(l)=FRf2):FRt2)=FR:GDT0522 

560 HV=NV(2):HP=HF-HV 

570 T4(1)=T3+500!T4(2)=T3+1000 

572 FDRI=IT02:T=(T3+T4(I))/2 

580 CP=. 24132+4. 4512E-54T“B.5504E~9IT[2+1.1B52E-12IT[3 

590 ER 1 1 )=HICP*IT4 ID-T31 -HP* 18200:NEXT 

600 IFABS (ER (2) ) < 1E-5THEH640 

610 T4=IER(lHT4(2)-ER(2)lT4(lll/IERtI)-ERI2)) 

615 T4(I]=T4(2):T4I2)=T4sBDT0572 

640 V4=N»RIT4/P4 

650 T5=T4I (V4/V5)HX-1I 

660 P5=P4ICV4/V5UX 

670 H2=P31 (V4-V3) 

680 N3=P4IV4I((V4/V5)I(X-I)-1)/(1-X) 

690 H4=P6t(V2-Vl) 

700 N5=PUIV1-V2) 

710 H= (HI +N2+N3+N4+lf5) IDF 
720 HEP=H/IIV1-V2)I144I 


ORIGINAL' PAGE IS' 
OF POOR QUALITY 


IHEP 1SFC IEFF I HP* 



ORtGtNAL PAGE (S 
730 HP«HlRPfl/660Q0 OF POOR QUALITY 

740 SFC«HFIRPHI30/I!P 
750 EFF-H/ (817.78) 

770 IPRINTUSING'IH III! 11. II Oil till ll.ll II, II 1.11 l.lll III. I l.llt 11,11 

H.II'jttH.RPH.fiF.CR.THiO, T5M40,PI/70.72il,P6/70.7264,NV/*iF,DF l «EP,SFC l EFf,HP 
780 101040 
790 END 

800 DATA 169,2404,284,69.6,48.4,26,32 
810 DATA 170,2600,284,69.6,48.6,26.30 
820 DATA 171,2600,252,69.6,48.6,27,74 
830 DATA 172,2600,252,49.4,48.6,27.82 
840 DATA 173,2596,225,49.6,48.6,29." 

850 DATA 174,2598,226,49.6,48.6,29.33 
060 DATA 175,2601,195,49.6,48.6,30.74 
870 DATA 174,2602,195,69.6,48.4,30.59 
800 DATA 177,2406,165,69.6,48.6,31.74 
890 DATA 178,2605,164,49.6,48.4,31.77 
900 DATA 179,2598,142,69,6,48.4,32.90 
910 DATA 100,2598,142,69.4,48,6,32.90 
920 DATA 190,2608,164,69.2,49.2,31,50 
930 DATA 191,2404,164,69.2,49.2,31.73 
940 DATA 192,2808,165,49.2,49.2,32.59 
950 DATA 193,2804,162,69.2,49.2,32,89 
960 DATA 194,2405,145,49.2,49.2,30.92 
970 DATA 195,2409,165,69.2,49.2,30.96 
980 DATA 196,2610,166,49.3,49.3,21.96 
990 DATA 197,2607,166,69.3,49.3,22.02 
1000 DATA 198,2809,166,69.3,49.3,24.80 
1010 DATA 199,2811,164,69.3,49.3,24.67 
1020 DATA 200,3005,145,49.3,49.3,25.37 
1030 DATA 201,3007,165,69.3,49.3,24.98 
1040 DATA 200,3005,165,49.4,49.4,25.39 
1050 DATA 201,3007,165,49.4,49.4,25.03 
1060 DATA 202,2810,165,69.4,49.4,25.39 
1070 DATA 203,2805,165,69.4,49.4,25.52 
1080 DATA 204,2599,146,49.4,49.4,25.45 
1090 DATA 205,2603,166,69.4,49.4,25.28 
1100 DATA 206,2601,146,79.4,47.9,27.94 
1110 DATA 207,2601,164,79.4,47.9,27,95 
1120 DATA 208,2602,166,69.4,44.4,23.72 
1130 DATA 209,2605,165,69.4,44.4,23.78 
1140 DATA 210,2599,146,59.4,44.4,19.15 
1150 DATA 211,2601,166,59.4,44.4,19.18 
1160 DATA 212,2405,165,80.5,49.5,26.71 
1170 DATA 213,2412,165,80.5,49.5,26.72 
1180 DATA 214,2605,165,67.5,44.5,23.27 
1190 DATA 215,2610,164,67.5,44.5,23.24 
1200 DATA 216,2598,195,69.0,44,5,22.85 
1210 DATA 217,2602,195,69.0,44.5,22.85 
1220 DATA 218,2600,224,72.0,44.5,23.19 
1230 DATA 219,2603,224,72.0,44.5,23.17 
1240 DATA 220,2600,255,75.2,44.5,23,35 
1250 DATA 221,2603,255,75,2,44.5,23,37 
1260 DATA 222,2602,286,78.4,44,5,23,43 
1270 DATA 223,2603,285,78.4,44.5,23.49 



DIESELA NOMENCLATURE 


ORIGINAL PAGE IS* 
OF POOR QUALITY 


AF A/F RATIO 

A1 A/F RATIO 

CP CONSTANT PRESSURE SPECIF HEAT 

CR COMPRESSION RA^IO 

CV CONSTANT VOLUME SPECIF HEAT 

DF DIAGRAM FACTOR 

EFF INDICATED EFFICIENCY 

ER ( I ) ERROR TERM 

El USED TO GET VOLUMETRIC EFFICIENCY 

FA F/A RATIO 

FRQ) FUEL FRACTION 

HP INDICATED HORSEPOWER 

M TOTAL MASS 

MA MASS DF AIR 

MEP INDICATED MEAN EFFECTIVE PRESSURE 

MF MASS OF FUEL 

MP MASS OF FUEL FOR CONSTANT PRESSURE PROCESS 

MV MASS OF FUEL FOR CONSTANT VOLUME PROCESS 

PI INTAKE PRESSURE 

PX PRESSURE TERM IN EXPONENT EQUATION 

P1-P6 PRESSURES 

Q TOTAL HEAT OF COMBUSTION 

QP Q USED DURING CONSTANT PRESSURE PROCESS 

QV Q USED DURING CONSTANT VOLUME PROCESS 

R GAS CONSTANT 

RN RUN NUMBER 

RPM REVOLUTIONS PER MINUTE 

SFC INDICATED SPECIFIC FUEL CONSUMPTION 

T1-T6 TEMPERATURES 

T4 ( I ) ESTIMATES FOR T4 

VE VOLUMETRIC EFFICIENCY 

V1-V5 VOLUMES 

W WORK SUM 

W1-W5 PROCESS WORK QUANTITIES 

X POLYTROPIC EXPONENT 

XI „ X2 TERMS IN POLYTROPIC EXPONENT EQUATION 



10 "DIESELS" 

20 'FUEL=C12H26 ORIGINAL PAGE IS! 

30 'FT=T0TAL HOLES FUEL:A=T0TAL UNITS AIR 0 p poOR QUALITY 

40 ’ F*=M0LES FUEL FOR CONST V0L:V=TDC VOL 
50 * T2=TEMP AT END OF C8HP:P3=1200PS1A 
60 DEFINT I,J,K,H,Z 

70 DIHAns t 16),B(15,16},XU6> l Y(i5) l SS{14,4),DO&),YRtl5) 

GO DIHL1(H),L2(14),NE(14J,HFU4),CV{141,URU41,HR(14) 

90 DIHHtl51,HlU4),U2(14),H3(14),H4ll4),H5U4),l)U5) 

100 DINC(14),Cl(14),C2(14) l C3(14) t C4(14} l C5{14) l SI(15) 

110 DIKStl4),Sltl4),S2U4),S3(14),S4U4),S5ll4},Q(9) f QQ(9) 

120 IN=16:N=15:X(16)=0:D!141=0!A=.01727sHF=”84.2 

121 LC=2.302585:R=1.98726:RB=. 082055 

125 INPUT“CR| T2-DES K, fl/F, P«AX-PSIA*j CR, T2, flF, Prt 

126 V=l/CR:P3=PH/14,696:FT=138IA/1170IAF) 

130 'READ PRODUCTS 

140 REST0RE:F0R3=1T015:READS413>:NEXT 

150 'READ NUMBER OF ATOMS OF C,H,M IN EACH SPECIE 

160 F0RJ=1TO14:FGR1=17Q4:READS8 (3,1): NEXTI ! NEXT3 

170 'READ COEFF FOR EQUIL EONS 

180 FOR I =1T09 J FOR J= 1 T01 4 : READA ( I , J > : NEXTO : NEXTI 

190 'READ CONST FOR EQUIL CONST EON 

200 F0RJ=1TO14:READL1(3>,L2(J):NEXT 

210 'READ THERHO-CHEN DATA: H=ENTHALPY, C=SPEC HEAT,S=ENTROPY 

220 F0RJ=1T014 

230 READNE(Jl,HF{J) l Hl(J},H2(J),H3(J) l H4fJ),H5(J} 

240 READCH3>,C2(3),C3(0),C4(J),C5I3) 

250 READS1 (J) , S2 IJ) ,83 ( 3) , S4 <J) ,S5 fJ) : NEXT 
260 'PROCESS: IF K=0 V=C0NST,K=1 P=CONST 
280 T3=1200:T4=2!00:F=1.22E-4 
290 K=0:T=T2:Z=0 

300 'INITIALIZE NO HOLES OF REACT FOR V=CUNST 
310 YRm=0:YR12)=0:YR13)=0:Yft(4)=3.76IA:YR(5)-0:¥R(6)=0 
320 YR{7>=0:YR(8>=0:YR<9)=A:YR(10>=0:YR!il>=0:YiM2>=0 
330 vn U3)=0:YR!14)=0 

340 „EAD NO HOLES OF PROD FOR FIRST ITERATION 
350 FOR J=1T01 4: READY (3) iNEXT 
360 'INTERNAL ENERGY OR ENTHALPY OF REACT 
370 UR=O:HR=0:YRU0)=O 
380 F0RJ=1T014 

390 HR(3)=100tNEl3)+HF(3)+Hl(3)+H2(3)IT+H313)ITI2 
400 HR (3)=HR ( J)+H4(JI ITE3+H5TJ) ITC4 
410 IFK=1THEN430 

420 URt3)=HR13)-RIT/1000:UR=UR+YR(3)IUR(3):G0T0440 
430 RR=HR+YR (J)tHRCJ) 

440 NEXT3 

450 HR=HR+F4(HF*2500) 

460 1 FK=0THENT=T3: Y 1 15 ) =T3: 60T0500: ELSET=T4 : Y ( 15 > =T4 
470 'READ NO HOLES OF PROD FOR FIRST ITER FOR P=CQNST 
480 F0RI=1T014:READY(I):NEXT 
490 ’TOTAL NO HOLES OF C,H,N,0 OF REACT 
500 YC=0:YH=0:YN=0:YO=0 

510 F0R3=1T01 4: YC=YC+SS (3 , 1 WR (3 ) : YH=YH*SS 13 , 2) IYR l 3 ) 

520 YN=YN+SS1J, 3) IYR (3 ) : Y0=Y8+SS ( 3 , 4 ) I YR ( J ) ; NEXT 

530 YC=YC+12*F: YH=YH+2MF 

540 ’TOTAL.NO HOLES OF C,H f N f O OF PROD 

550 C=0: H=0: N=0 : 0=0 : Y ( 1 0) =0 : YP=0 

560 F0RJ=1T014 

570 C=C+SS (3,11 *Y 13) :H=H+SS (3 ,2) tY (3> : K=N*-SS 13,3) *Y 13 ) 

580 Q=OtSS 13, 4)IYI J) :YP=YP+Y 1 J) :NEXT 



590 ’PRESSURE 

400 P=YPtR8IT/V 

610 ’CQEF OF ATOM BAL EONS 

620 F0RI=10TO13.‘FORJ=lT014:A(I,J)=SS(J,I-9)JY(J)iNEXTJ:NEXTI 
630 ’ COEF OF XF FOR ATOM BAL AND PRESS EONS 
640 1FK=1THEN670 

650 A(10,10)=-12IFIC/YC:AU1,10)=-26IFIH/YH 

460 A(12,10)=0:A(13,10)=0:A(14,10)=0:60T0690 

670 AUO,10! s 0:A(ll,10)=0:A(12,10)=OiA(13,10)=0iA(l5,10)=Q 

6S0 ’COEF OF PRESS EON 

690 F0R3=lT014:A(14 f J)=YU):NEXT 

700 'ERROR CONSTANTS 

710 D ( 10) s L0S (C/YC) /LC; D ( 1 1 ) =LOS IH/VH) 7LC 

720 D(12)=LOStN/YN)/LC:D(13)=L0S(0/Y0)/LC:Dll4)=L0B(P/P3)/LC 

730 A(iO,16)=-CIDtiQ)sA(ll,16)=-HID(ll) 

740 A(i2 t 16)=-NID«2J;A113,16»=-0ID113):AU4,16)=“YPID(14) 
750 NP=0:UP=0:SP=0 

760 ’ENTHALPY OR INTERNAL ENERGY OF PROD 
770 F0RJ=IT014 

7B0 H(J)=100*NE4J>*HF(J>+H1 ( J ) +HZ4J) IT+H34J) ITt2 
790 H4J)=HtJ)+H4(J)*TE3+H5(J>ITE4 
800 G0T0840 
BIO ’ENTROPY 

820 S(J)=Sl(J)+S2(J)mS3!J)m2+S4(J)ITE3+S5{J)4T[4 
030 SiJ)==(S{J)-RtLOG(Y(J)*RBIT/VI)/1000 
840 IFK=1THEN870 

B50 U(J)=H(J)-RIT/1000:UP=UP+YIJ!IU(J) 

860 G0T0889 
870 HP=HPtY(J)IH(J) 

880 HEXTJ 

890 ’COEF OF ENERGY ESN 
900 IFK=1THEN920 

910 F0RJ=lTD14r A (15 f J)-Y(J)IU0J) sNEXTi G0T0960 
920 F0R3=lT014:A(15 f J)=Y(3)IH(J):NEXT 
930 'COEF OF XV IN PRESS EON FDR P=CONST 
940 AC14,10)=-YP 

950 ’COEF. OF XT IN PRESSURE EQN. 

960 A(14,15)=YP 

970 'INTERNAL ENERGY OF PRODUCTS FOR V=C0NST 

980 IFK=1THEN1030 

990 UF=NF+2500: i)P=UP-F4UF 

1000 ’COEF OF XF IN ENERGY EQN FOR V=CQNST 

1010 A(15,10)=-UF1F:GQTD105Q 

1020 'COEF OF XV IN ENERGY EQN FOR P=C0NST 

1030 A (15, 10) =0 

1040 ’COEF OF XT AND XV IN EQUIL EQN 
1050 QD=RIT/10Q0 

1060 FDRI=1T09;A(I,15)=0:A(I,10)=0 

1070 FORJ=1T014 

1080 ’COEF OF XT IN EQUIL EQN 

1090 A(I,15)=AU,15)+A(I,3)-A(l,J)lHm/QD 

1100 1FK=0THEK1 I 20 

1110 A(I, 10)=A(I,10)-A(I,J) 

1120 NEXTJiNEXTI 
1130 CV=0:CP=0 

1140 ’COEF OF XT IN ENERGY EQN 
1150 FDRJ=1TQ14 

1160 C(J)=Cl(J)iC2(J)mC3(J)m2tC4(J)IT[3+C5(J)*T[4 
1170 IFK=tTHEN1190 

1 180 CV ( J) =C ( J ) -R: CV=CV+Y ( J) ICV t J ! ! B0T01200 
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1190 CP*CP+YCJ)IC(J) 

1200 NEXTJ 

1210 IFK*1THEN1240 ORIGINAL PAGt; f;I 

1220 0 ( 15, 1 5) =CVIT/ 1000 1 D ( IS) “LOG (UP/UR) /LC OF POOR QUALITY 

1230 A (15 t l6)«-*llPtD(15) i BDTD1270 

1240 Al 15, 15)=CPH/ 1000: DUS) =lOli (HP/HR) /LC 

1250 A(15,t6)*-HPiDU5> 

1260 ’ERROR CONSTANT FOR EQUIL EON 
1270 F0RI=1TQ9 

1280 Btn=LOG(YU>)-SStI,niLQG(Y(14))-SS(I l 2)tLOS(YtlJ>) 

1290 B(I}=B(I)-SS(I,3HLQG(Y(12))-SSU,4)ILGGIYU3H 

1300 fl(n=B(I)Ml-SStI,l)-SSU,2)-SSII,3l*5SII,4)) IL0GIRBIT/V) 

1310 B0(l) s LltI)“L2II)/T 

1320 flO(I)=fiO(II-S5(],llt(LH14)-L2U4)/T)~SS(I,2) I tLM 1 D-L2U 1) /T> 

1330 B0(I)=B0(I)“SSU,3>IIL1(12I-L2U2)/T)-SSII,4) I(L1I13)-L2U3)/TI 
1340 DtI) s GM) /LC-Q04I } : A( I* 16) s -0 (I) (NEXT 
1350 EP=Q:EB=0 

1360 FORI=1TO 3:EC=EQ+ABS(D(I))iNEXT 

1370 FORlMTOlSiERsEP+ABStDinhNEYT 

1380 IFK=0THEN1400 

1390 EP=EP+EQ 

1400 PRINTEP,T 

1410 IFEP<. 0005THEN1750 

1420 AU0,15)=0!A(U,15)=0;AU2,15)s0!Atl3,15)=0 
1430 'MATRIX SOLUTION 

1440 FORI=lTQMs XIII =0:FDRJ=1TOIN:B ( I, J) =0:NEXTJs NEXTI 

1450 FQRI=1T0H:B(I,1)=AM,1)!NEXT 

1460 F0RJ=2TQINiB(l,J)*AU,J)/Atl,lhNEXT 

1470 F0RIC=2TQK:J=IC 

1480 FGRl»JTOH:SUH=0:Jl=iM 

1490 FQRIK=lTDJlsSUH=SUH+B(I, IK) #B ( IK, sNEXTIK 

1500 BU,J)sAII,J)-SUflsNEXTI 

1510 I=J:I1=IU 

1520 FORJ=IlTOINjSUN)iJl=J-i 

1530 FQRlK=tTQJliSUH*SlllRBU,IK)tB(IK,J):NEXTIK 

1540 B(I,J}=IA(I,J)‘SU|1)/B{I,!>:NEXTJ:NEXTIC 

1550 IFZ<3TRENL=.06ELSEL=. 1 

1560 X tM)=B(M, tM) t IFX4M) >LLETX (M)=L 

1570 IFX(HK-LLETX(ft)=-L 

1580 Y(H)=Y (M) tlOEX 4M) : IFY IN) <900LETY M) =900 

1590 IFY U1) >3500LETY (N) =3500 

1600 MIN) 

1610 FORIL*2TQN: M-IL+N: SUM=0:I1=I+1 

1620 FORIK-I 1T0N: SUM=SUM+B II , IK) IX t IK) t NEXTIK 

1630 X 1 1 ) =B 1 1 , 1 N) -SUN: IFX II ) 736LETX 1 1 ) =36 

1635 IF X(IK-36LETXII)=-36 

1640 YII)=YIIU10[XU):IFYUKlE"38LETYII)=lE-38 

1650 IFY(I)>3.76IALETYII)=3.76IA 

1660 NEXTIL 

1670 Z=Z+1:IFK=1THEN1690 
1680 F=FI10IX (10) : IFF<FT/50LETF=FT/50 
1685 IFF>FTLETF=FT 
1687 60T0550 

1690 V=VI101XU0):IFV<1/CRLETV=1/CR 
1700 IFV>1LETV=1 
1710 G0T0550 
1750 YP=O:XP=0 

1760 FOR J=1T0X4! YP=YP+Y { J ) : NEXT 

1770 F0RJ=lT014:X(J]=YIJ)/YP:XP=XPU(J)iNEXT 

1780 IFK=1THEN1970 



1785 lPRIKT , CR«*jCR, , T2= , ;T2, , fi/F=*jftF,'PHft7* , |PH 
1790 LPRINT'CQNSTANT VOLUME PROCESS' 

1800 LPRINT'INPUTjNO HOLES OF FUEl a '|FT 
1810 LPR1NT TAB (61 'HO HOLES OF AIR*'}AI4.76 
1820 LPRINT TAB (6) 'PEAK PRESS P3 »'jP|*ATn* 

1830 LPRINT TAB ti) •TEMPERATURE T2=';T2 
1810 LPRINT*V=NO OF HOLES OF PRODUCTS* 

1050 LPRINT* X=HOL£ FRACTION OF PROOUTS' 

I860 FOR 0=1TQ14 
1870 LPRINT *Y("jS4(J);*)='jY(J)}TAB(21)'Xl*;Sf (J) j*) =, jXlJ) 

1880 HEKTJ 

1890 LPRINT’YT0T c "iYP}TAB(211 , V0LUHE= , , , V 

1900 LPRINT'TEHP T3=' | T? TAB (21 ) ■UR='jUR 

1910 LPRINTTJEL=*jFjTAB(21)*UP='}UP 

1920 LPRINT' PRESSURE 5 ' } P } T AB (21 ) ' 1 TER AT ION= ' j 2 

1930 LPRINT'TOTAL ERADR=" f EP j LPRINT* '.'LPRINT* ' 

1940 ’INITIALIZE NO HOLE OF REACTANT FOR P=CONST PROCESS 
1950 F0RJ«1T0I4i YR(J) b Y (J) .'NEXT 
1960 F=FT-F;K=liZ=0:V=.18:GGTQ370 
1970 LPRINT'CONSTANT PRESSURE PROCESS' 

1980 F0R3=1T014 

1990 LPRINT'Yt'jS4(i]}j'} :: "jY(J)(TAB(2irX('}St{Jl;*)=*{X(J} 

2000 NEXT! 

2010 LPRINT*YTQT=';YP;TAB(21)'VQLUHE=';V 
2020 LPRINT'TEHP T4=' j T? TAB C21 > *HR=' ;HR 
2030 LPRINT'FUEL=*;FjTAB(21)*HP='|HP 
2010 LPRINT" PRESSURE =';P;TAB121} , ITERATI0N=*;Z 
2050 LPRINT'TOTAL ERRORS }EP 
2055 60T0125 

2070 DATA CQ2,C0,CN4,N2,N0,H2Q,H2,0H,02,F,H,N,Q,C,T 

2080 DATA 1,0, 0,2,1, 0,0, 1,1, 4,0,0, 0,0, 2,0, 0,0, 1,1, 0,2, 0,1 

2090 DATA 0,2, 0,0, 0,1, 0,1, 0,0, 0,2, 0,0, 0,0, 0,1, 0,0, 0,0, 1,0 

2100 DATA 0,0, 0,1, 1,0, 0,0 

2110 DATA 1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,-2, -1 

2120 DATA 0,1, 0,0, 0,0,0, 0,0, 0,0,0, -1,-1 

2130 DATA 0,0, 1, 0,0, 0,0, 0,0,0, -4,0, 0,-1 

2140 DATA 0,0, 0,1, 0,0, 0,0, 0,0,0, -2, 0,0 

2150 DATA 0,0,0,0,l,0,0,0,0,0,0,-i,-l,0 

2160 DATA 0,0, 0,0, 0,1, 0,0, 0,0, -2,0, -1,0 

2170 DATA 0,0, 0,0,0, 0,1,0, 0,0, -2, 0,0,0 

2180 DATA 0,0, 0,0, 0,0, 0,1, 0,0, -1,0, -1,0 

2190 DATA 0,0,0, 0,0, 0,0, 0,1, 0,0,0, -2,0 

2200 DATA -.0295997,-20764.8,4.2918,-6345.6,-5.8028,-4790.4 

2210 DATA 0, 0,. 6574, 4711.2,-3.0466, -13168.8,0,0,. 7068, I8B6.4 

2220 DATA 0,0,0,0,3.1626,11896.8,3.514,25116 

2230 DATA 3.503,13356,8.1622,37281.6 

2240 DATA 2, -94.054, -4.04934,. 0106409, 1.66486E-6, -2.97619E-10,2.1205SE-14 

2250 DATA 8.65961, 6.62203E-3r2.80545E-6,5.50571E-10, -4.04973E-14 

2260 DATA 47. 3063,. 0224858, -6.43573E-6, 1.06715E-9, -7.18642E-14 

2270 DATA 1.5, -26.417, -2.21792 I 6.65945E-3,B.B95E-7, -1.60617E-10,1. 14865E-14 

2280 DATA 5.81533, 3.20923E-3,-1.32B46E-6,2.55813E-10, -1.85551E-14 

2290 DATA 45.6148, .0138031, -4. 00914E-6, 6. 70106E-1Q, -4.53569E-14 

2300 DATA 3,-17.895,-4.128l9,7.B5264E-3,6.44417E-6, -1.13465E-9,7.77547E-14 

2310 DATA 2 . 29288 , , 0222002, -8 . 84 428E-6 , 1 . 64 1 38E-9 , -1.16043E-13 

2320 DATA 37.503B,.0268B12,-6.04799E-6 I 8.54565E-10, -5.18365E-14 

2330 DATA 1,0, -2.08601, 6. 39974E-3,9,78717E-7,-1.7702E-IO, 1.26646E-14 

2140 DATA 5.61035, 3.31765E-3,-1.34BB3E-6,2.56255E-10, -1.B4078E-14 

2350 DATA 44. 0463,. 0139597, -4.16865E-6,7.I8432E-10, -5.00693E-14 

2360 DATA t, 21.58,-2.3971, 7.10077E-3.7.27649E-7, -1.29337E-10,9.13029E-15 

2370 DATA 6.2401, 2.87745E-3,-1.20965E-6,2.35588E-I0, -1.72297E-14 
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2JB0 DATA 48.584,. 0144177, -4. 30522E-6,7.35191E-10, ~5.O6504E-14 

2390 DATA l.5,-S7.79B,-1.74099,5.6689BE-3,2.4t?4E-6, -4.16138E-10,2.65073E-14 

2400 DATA 5.21619, 6.19274E-3,-1.79861E-6,2.5224lE-10, -1.36954E-14 

2410 DATA 43.2465, .0156744, -3. 89475E-6, 6. Q8863E-10, -4.00555E-14 

2420 DATA 1,0, -1.39667, 5.4312E-3,l.04227E-6,-1.4l029E-l0, 9.13029E-15 

2430 DATA 4.00244, 1.118fl9E-3,7v44342E-a, -7. 15413E-11, 8.32035E-15 

2440 DATA 30.5792,. 0119128, -3,2945E-6,5.48002E-10, -3.71102E-14 

2450 DATA I, 9.33, -1. 57142, 5.49288E-3,1.01526E-4, -1.44545E-10.8.B3577E-15 

2440 BATA 5.72447, 1.99422E-3,-4.27002E-7, 3, 39192E-11,0 

2470 DATA 43.0974,. 0123579, -3.47U7E-4,5.79225E-10, -3.94444E-14 

2480 DATA 1,0, -2.40987, 7.44459E-3,3.84913E-7,-1.21984E-ll, -8.83577E-16 

2490 DATA 7.12922, 1.42729E-3, -5. 17909E-7,1.09944E-10, -9.5720BE-15 

2500 DATA 47.1104, .0148438, -4.51712E-4.7.94175E-10, -5.5959BE-14 

2510 DATA 0,0, 0,0, 0,0,0 

2520 DATA 0,0, 0,0,0 

2530 DATA 0,0, 0,0,0 

2540 DATA .5, 52.102, -1.48109, 4.94807E-3,-1.21905E-ll, 0,0 
2550 DATA 4.948,0,0,0,0 

2540 DATA 24.6523, 9.23922E-3,-2.9799BE-4,5.3146E-10, -3.74993E-14 

2570 DATA .5, 112.945,-1. 45182, 4.B905SE-3,7.43347E-8, -3.0B721E-11,4.71241E-15 

2580 DATA 4.84307, 2.16256E-4,-1.453Q3E-7,3.50204E-il, -1.693S2E-I5 

2590 DATA 35.95B3,9.0aOA2E-3,-2,B7709E-6,5.03335E-10, -3.4754E-14 

2400 DATA .5,59.559,-1. 47492, 5.06465E-3.-4.23555E-8, 4.02427E-12,2.94524E-16 

2410 DATA 5,03035, -1 . 64497E-5, -2. 8 1 153E-8, 1 . 40049E-1 1 , -1.25173E-15 

2420 DATA 37.9917, 8. 99493E-3,-2.81594E-4,4.88005E-10, -3.35759E-14 

2630 DATA 1,170.886, -1.52089, 5.11104E-3,-1.32952E-7, 4.5I027E-U,-3.38704E-15 

2640 DATA 5.09745, -2.2207E-4,9.69728E-8,-i.l5443E-12 r -1.32537E-15 

2650 DATA 37.1575, 9.Q0799E-3,-2.B5048E-6,5.04278E-lQ, -3.53431E-14 

2640 DATA 1.5E-3,7.5E-12, IE-38,. 065, 1.7E-5,i,6E-3,4.BE-12 

2670 DATA 3.6E-8, . 015, 0, 1 . 2E-14 , 3.0E-20, 9. 0E-1 1 , IE-38 

2680 DATA 6.4E-3,7.3E-6,1.0E-22,.065,5.5E-4,7.1E-3,1.7E-6 

2690 DATA 6.9E-5,6.9E-3,0,1.3E-7,3.0E-lt,3.7E-6,3.6E-23 

2700 END 



DIESELB NOMENCLATURE 


A 

A ( I , J ) 

AF 

B ( I , J ) 

C 

C(J) 

CP 

CR 

CV 

CV(J) 

C 1 *-C5 ( J) 
D ( I ) 

EP 

EG 

F 

FT 

H 

H ( J) 

HF 

HF ( J ) 

HP 

HR 

HR ( J ) 
H1-H5(J) 
I, IC, IK 
III, IN, II 

j*,n,k 

L 

LC 

LI -L2CI) 

M 

N 

NE ( J) 

0 

P , P3 
Q ( I ) 

QD 

GO < I ) 

R, RB 
S(J) 

SP 

SS(I,J) 
S1-S5 (J) 
S$ (J) 

T, T2-T4 
U(J> 

UF 

UP 

UR 

URCJ) 

V 

X a> 

XP 

Y ( J ) 

YC, YH 
YN, YD 


TOTAL UNITS OF AIR (02+3.76 N2> 

MATRIX ELEMENTS 
AIR FUEL RATIO 
MATRIX ELEMENTS 
C-ATOMS IN PRODUCTS 
SPECIFIC HEAT OF PRODUCT J 
SPECIFIC HEAT OF PRODUCTS 
COMPRESSION RATIO 
SPECIFIC HEAT OF PRODUCTS 
SPECIF HEAT OF SPECIE J 
SPEC HEAT EQUATION CONSTANTS 
DIFFERENCE FACTOR 
ERROR PARAMETER 
ERROR PARAMETER 

MOLES FUEL FOR CONSTANT VOLUME PROCESS 

TOTAL MOLES OF FUEL 

H-ATOMS IN PRODUCTS 

ENTHALPY OF PRODUCT J 

ENTHALPY OF COMBUSTION 

ENTHALPY OF FORMATION OF SPECIE J 

ENTHALPY OF PRODUCTS 

ENTHALPY OF REACTANTS 

ENTHALPY OF REACTANT J 

ENTHALPY EQUATION CONSTANTS 

INTEGERS 

INTEGERS 

INTEGERS 

ERROR LIMIT 

LOG CONVERSION CONSTANT 
EQUILIBRIUM EQUATION CONSTANTS 
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INTEGER 

N-ATOMS IN PRODUCTS 

MOLES OF ELEMENTS IN SPECIE J 

Q-ATDMS IN PRODUCTS 

PRESSURE 

CALCULATED EQUIL CONSTANT FOR FORMATION OF I 
RJkT/lOOO 


EQUILIBRIUM CONSTANT FROM JANAF TABLES 
GAS CONSTANTS 
ENTROPY OF PRODUCT J 
ENTROPY OF PRODUCTS 

NUMBER OF ATOMS OF TYPE I IN SPECIE J 
ENTROPY EQUATION CONSTANTS 
SPECIE SYMBOL 
TEMPERATURES 

INTERNAL ENERGY OF SPECIE J 

INTERNAL ENERGY OF FUEL 

INTERNAL ENERGY OF PRODUCTS 

INTERNAL ENERGY OF REACTANTS 

INTERNAL ENERGY OF REACTANT J 

TOP DEAD CENTER VOLUME 

MATRIX SOLUTIONS 

TOTAL MOLE FRACTION OF PRODUCTS 

MOLES OF PRODUCT J 

C- AND H- ATOMS IN REACTANTS 

N- AND 0- ATOMS IN REACTANTS 



VP 

YR< I) 
Z 


MOLES 2° L Rl fl ?U? D ? UCTS 
number of iterations 
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A/F* 24 


PMAX* 1200 


CR= 10 T2* 777 

CONSTANT VOLUME PROCESS 
INPUT : NO MOLES OF FUEL' 5.84132E-04 
NO MOLES OF AIR= .0022052 
PEAK PRESS P3 = 81.6549 ATM 
TEMPERATURE T2= 777 
Y=NO OF MOLES OF PRODUCTS 
X=MOLE FRACTION OF PRODUTS 


Y <C02)= 1.46337E-03 
Y (CO) “ 7. 51498E—12 
Y (CH4) = IE— 38 
Y(N2)= .0649267 

Y ( NO ) * 1 . 66875E-05 
Y (H20) = .0015853 
Y ( H2 ) = 4. 77274E-12 
Y (OH) = 3. 63892E-0B 
Y (02) = .0150056 
Y(F>= 0 

Y <H) - 1.21079E— 14 
Y(N) = 2. 99162E— 20 

Y <0> = 9. 01939E— 1 1 
Y(C> = 1E-3B 
YTOT* ,0829977 
TEMP T3= 1198.98 
FUEL* 1.21754E-04 


X (C02> = .0176314 
X(CO>* 9. 05444E— l 1 
X (CH4)= 1.204Q5E-37 
X <N2> = .782271 
X (NO) = 2. 0A06E-04 
X (H20) = .0191005 
X (H2) = 5. 75045E— 1 1 
X (OH) = 4. 38437E— 07 
X (02) — .180796 
X CF) = 0 

X CH> = 1 . 45B83E-13 
X <N> = 3. 60446E— 19 
X<0) = 1 . 0B67E-09 
X (C) = 1 . 20485E— 37 
VOLUME* . 1 
UR* 8.37599 
UP* 8.37599 


PRESSURE* 81.6549 ITERATION* 3 
TOTAL ERROR* 4. 134SE-04 


CONSTANT PRESSURE PROCESS 
Y (CG2) * 6. 99B39E-03 
Y (CO) * 1.35258E-05 


Y (CH4) = 6. 65312E-22 
Y(N2)= .0646316 
Y (NO) * 6.06853E-04 
Y ( H20 ) = 7. 54721E-03 
Y ( H2 ) * 2. 92444E-06 
Y (OH) = 9. 1 9496E—05 
Y ( 02 ) * 6. 13907E-03 

Y < F ) « 0 

Y (H) - 2.61932E-07 

Y (N) * 6. 9645 IE-1 1 
Y (0) = 5. 57688E-06 
Y(C)= 3.09552E-22 
YTOT* .0860374 
TEMP T4= 2167.47 


X (C02) * . 0813413 
X(CO)= 1 . 572O0E— 04 
X (CH4) = 7. 73203E-21 
X (N2) = .751204 
X <N0) * 7 . 05336E— 03 
X (H20) * .0877201 
X ( H2 ) = 3. 39903E— 05 
X <0H) * 1 . 06872E— 03 
X (02)= .0713535 
X (F) = 0 

X (H) = 3. 04439E— 06 
X (N) * 8. 09476E—10 
6. 4S192E—05 
3. 59788E— 21 
. 187398 


X CO)* 

X <c> = 
VOLUME* 

HR* 9.9849 
FUEL* 4. 62379E— 04 HP* 9.98489 

PRESSURE * 81.6549 ITERATION* 
TOTAL ERROR* 3. 19628E-05 
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